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摘要 :黄土 高 原 退 耕 还 林 ( 草 ) 等 生态 恢复 工程 促进 了 地 表 植 被 覆盖 增加 ,进而 通过 影响 地 表 - 太 气 之 间 热 量 交换 影响 区 域 气候 
过 程 。 基 于 黄土 高 原 1998 一 2000 年 和 2008 一 2010 年 SPOT 卫星 反 演 的 植被 覆盖 资料 ,54 个 地 面 气 象 站 气温 资料 及 EIN-Interim 
地 表 热 通 量 数据 ,采用 空间 分 析 交 叉 验 证 及 地 表 热 量 平 衡 分 析 的 方法 ,从 站 点 尺度 探讨 了 退耕 还 林 ( 草 ) 工 程 初期 和 10 年 后 黄 
土 高 原 植 被 变化 与 气温 和 地 表 热 通 量变 化 之 间 的 关系 。 研 究 表明 ,退耕 还 林 ( 草 ) 工程 开展 10. 年 后 ,黄土 高 原 气 温 最 小 值 、 最 
大 值 与 平均 值 均 有 下 降 , 植 被 覆盖 增加 与 气温 变量 降低 在 空间 上 呈正 相关 。 同 时 ;植被 覆盖 增加 与 潜 热 通 量 增加 `, 感 热 通 量 与 
大 气 下 行 长 波 辐射 下 降 在 空间 上 也 呈正 相关 关系 。 这 些 结果 表 明 ,植被 恢复 可 通过 增加 地 表 蒸 散发 作用 对 区 域 气候 产生 降温 
效应 ,会 减缓 气温 升 高 对 黄土 高 原生 态 系统 的 影响 。 
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Abstract: A large-scale re-vegetation with regional ecological engineering on the Loess Plateau, including the Grain for 
Green Project (GGP), has already greatly impacted the heating and energy flux between the land and atmosphere, which 
induces regional climate. We used two methods, the cross-checking of spatial change and surface heat balance analysis to 
detect the relationships aniong the vegetation change, temperature change, and surface heat flux change during the GGP at 
the station scale over the Loéss Plateau. The 10 days normalized difference vegetation index ( NDVI) derived from SPOT 
data was used’ to indicate-the vegetation. The monthly temperature data from $54 weather stations within or nearby the Loess 
Plateau weére~used to describe the temperature change. The EIN-Interim surface heat fluxes were analyzed to determine the 
surface heat condition. Two stages, the initial stage (1998—2000, Stage IT) of GGP and 10 years after ( 2008—2010, 
Stage I[) were compared. The minimum, maximum, and mean temperatures were lower in Stage Il, and the increase of 
vegetation was spatially positively correlated to the decrease in all temperature variables. Meanwhile, the increase in 
vegetation was spatially positively correlated to the increase in latent heat flux, and the decrease of sensible heat flux and 
the atmospheric downward longwave heat flux. These results indicated that vegetation restoration may lower the temperatures 


through an increase in surface evapotranspiration, and thereby mitigate the climate warming effect on the ecological system 
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of the Loess Plateau. 
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植被 作为 陆地 生态 系统 的 主体 ,也 是 气候 系统 中 重要 的 动态 的 组 成 部 分 。 陆 地 植被 可 以 通过 改变 
地 表 反 照 率 及 陆地 表面 净 辐 射 在 感 热 潜 热 与 地 面 热 通 量 之 间 的 分 配 “ ;或 通过 改变 表面 粗糙 度 “ 来 调节 大 
气动 力 特 征 , 进 而 调 市 地 - 气 间 的 水 分 和 能 量 循环 ,影响 地 表面 能 量 平衡 ,从 而 影响 气候 特征 。 目 前 在 局 地 和 
区 域 的 尺度 探寻 植被 对 气候 的 影响 过 程 与 机 制 方面 已 经 开展 了 众多 人 研究“ 。 但 是 ,由 于 地 球 表 面 物 理化 
学 和 生物 特性 存在 明显 区 域 差异 ,植被 对 气候 要 系 影响 的 机 制 也 可 能 会 因 区 域 不 同 而 存在 差 寞 。 例 如 , 侈 伐 
森林 可 能 导致 热带 地 区 的 气温 升 高 ,但 在 寒冷 地 区 气温 反而 下 降 …”。 因 此 ,为 了 进一步 认识 植被 与 气候 因子 
之 间 的 关系 ,有 必要 再 在 不 同 地 区 和 不 同时 空 太 度 下 分 析 植 被 分 布 和 变化 过 程 .地 表 参 数 及 气候 对 其 的 响应 。 

页 土 高 原 是 中 国 气 候 变化 的 敏感 区 域 之 一 ,其 独特 的 目 然 气候 状况 以 及 生态 环境 变化 受到 中 外 学 者 广泛 
关注 “”。 为 缓解 中 国 北方 地 区 的 水 土 流失 和 土地 退化 等 问题 , 自 1999 年 开始 在 黄土 高 原 开展 了 大 规模 的 退 
耕 还 林 ( 草 ) 等 生态 工程 ,植被 增长 效果 显著 。 目 前 对 黄土 高 原 退 耕 还 林 ( 草 ) 政 策 下 的 农户 特征 变化 、 
经 济 效益 ”研究 较 多 ,而 对 于 该 政策 区 域 环境 效应 ,尤其 是 植被 恢复 后 其 与 大 气 之 间 的 相互 作用 的 探讨 较 
少 。 郭 建 侠 等 使 用 中 尺度 气象 模式 耦合 陆 面 植被 模式 模拟 退耕 还 林 对 局 地 环境 要 素 的 影响 ” ; 王 平等 基于 
CLM 模型 对 黄土 高 原 植 被 覆盖 变化 对 气温 和 降水 的 影响 作 了 模拟 研究 ” ; 梁 玲 等 使 用 中 尺度 气候 模式 模拟 
了 黄土 高 原 植被 改变 对 夏季 气候 影响 “ ; 张 静 辉 等 使 用 区 域 气候 模式 研究 黄河 源 区 植被 覆盖 度 对 区 域 气候 
影响 “ 。 但 以 上 研究 大 多 基于 模型 模拟 的 方法 , 较 少 使 用 观测 数据 证 明 局 地 或 区 域 尺度 上 植被 恢复 对 气候 
的 反馈 作用 。 对 于 黄土 高 原来 说 ,基于 观测 数据 的 退耕 还 林 ( 章 ) 工程 是 否 可 能 导致 区 域 气温 变化 尚未 清楚 。 
本 文 利 用 遥感 卫星 反 演 的 植被 覆盖 数据 ,结合 54. 个 位 于 黄土 高 原 还 林 ( 草 ) 工 程 区 的 气象 站 资料 及 EIN- 
Interim 热 通 量 数 据 , 采 用 空间 分 析 交 又 验 证 及 地 表 热 量 平衡 分 析 的 方法 ,从 站 点 尺度 探寻 三 者 间 的 可 能 关 
系 , 分 析 植 被 变化 对 气温 影响 的 基本 过 程 与 其 物理 机 制 。 


1 研究 区 与 数据 


黄土 高 原 地 区 位 于 黄河 中 上 游 地 区 , 占 地 面积 大 约 62.4x10” km?, 介 于 33°43’ 一 41°16'N,100°54’ 一 114° 
33 守 之 间 , 包 括 太行 山 以 西 ,青海 日 月 山 以 东 , 秦 岭 以 北 阴山 以 南 广大 地 区 "(图 1)。 该 区 土壤 深厚 ,主要 土 
壤 类 型 为 褐 土 . 黑 坊 土 和 黄 绵 土 ,这 些 土壤 毛管 孔 过 大 ,易于 芋 发 ,保水 性 较 差 ." ;而 且 该 区 为 典型 大 陆 性 季 
风气 候 , 大 部 分 地 区 平均 年 降水 量 400 mm 左右 ,空间 分 异 较 大 ,而 且 季 节 性 分 布 明 显 , 约 有 69.33% 降 雨量 以 
暴雨 形式 集中 在 ,6 一 9 月 。 黄 土 高 原 的 独特 地 形 地 貌 结 构 及 其 土质 玻 松 和 暴 雨 频 发 等 特征 ,导致 该 地 水 土 流 
失 十 分 严重 ,是 我 国生 态 脆弱 区 和 水 土 保 持 重点 区 。 

本 文 主要 采用 基于 地 面 气象 站 的 气温 观测 数据 .基于 遥感 卫星 反 演 的 植被 NDVI 数据 与 EIN-Interim 再 分 
析 资 料 3 类 研究 数据 。 所 用 地 面 实际 观测 资料 来 自 中 国 气象 数据 共享 网 1998 一 2000 年 .2008 一 2010 年 黄土 
高 原 站 点 的 气候 资料 月 值 数据 集 。 

在 黄土 高 原 及 周边 地 区 选取 了 54 个 气象 站 点 ,选取 原则 为 :(1) 站 点 气温 数据 连续 ,在 1998 一 2000 年 、 
2008 一 2010 年 5 月 至 9 月 无 缺失 情况 ;(2) 站 点 位 于 还 林 还 草 工程 的 覆盖 区 ,保证 站 点 周 于 在 退耕 还 林 工 程 
后 期 主要 土地 类 型 为 草地 或 林地 (绿地 占 比 大 于 33.33%)。 选 取 的 站 点 分 布 如 图 1 所 示 , 气 象 站 10 km 绿地 
比重 是 根据 黄土 高 原 2005 年 1:25 万 土地 利用 数据 在 ArcGIS 空间 分 析 平 台中 计算 得 到 。 

NDVI 数据 集 由 SPOT 卫星 搭载 的 植被 传感器 SPOT VGT 反 演 得 到 ,空间 分 辩 率 为 1 km, 使 用 每 月 10 d 
最 大 值 合 成 的 NDVI 产品 (VGT-S10)。1998 一 2008 年 VGT-S10 数据 来 源 于 “黑河 计划 数据 管理 中 心 ” 
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图 1 研究 区 位 置 与 气象 站 分 布 


Fig.l Location of the study site and .distribution of weather stations 


( http://westdc.westgis.ac.cn) ;2008 一 2010 年 使 用 比利时 VITO 研究 所 图 像 处理 与 存档 中 心 提供 的 VGT-S10 
数据 产品 (http://www.vito-eodata.be) 。 本 文选 取 1998- 一 2000 年 ,2008 一 2010 年 ERA-Interim 月 平均 地 表 热 
通 量 资料 。ERA-Interim 数据 集 的 水 平分 辨 率 为 0.1T255x0.125° ,垂直 分 辨 率 60 层 ,ERA-Interim 再 分 析 资 料 同 
化 方案 .数据 源 等 参照 相关 文献 “ 


2 方法 


2.1 数据 处 理 

气温 数据 为 植被 生长 季 (5 一 9 凡 )54 个 气象 站 的 月 平均 气温 (了 ,Mean temperature) 月 最 高 气温 (了 ，， 
Maximum temperature) 与 月 最 低 气 温 ( 7 ，,Minimum temperature ) 的 均值 。NDVI 月 值 为 基于 VGT-S10 产品 的 
DN 值 转换 的 .NDVI 数据 的 月 最 大 值 合成 值 。 选 取 ERA-Interim 数据 集 的 植被 生长 季 感 热 通 量 (SH ,Sensible 
heat flux) 、 潜 热 通 量 ( LH ,Latent heat flux) 与 大 气 下 行 长 波 辐射 (LD ,Atmospheric downward longwave radiation ) 
3 个 热 通 量 值 分 析 地 表 热 量 平衡 状况 。 计 算 生 长 季 热 通 量 均值 并 使 用 克 里 金 插值 的 方法 估计 54 个 气象 站 的 
辐射 量 值 汪 所 有 的 站 点 气温 值 ` NDVI 值 及 热 通 量 采 取 站 点 周围 9 个 像 元 平均 值 。 
2.2\ (气温 、NDVI 的 时 空 变化 特征 

对 1998 一 2000 年 和 2008 一 2010 年 两 个 时 段 的 植被 生长 季 NDVI( NDVI in growing season， 以 下 记 为 
NDVIg) 和 气温 变化 的 空间 格局 进行 了 分 析 。 同 时 利用 AreGIS 10.0 的 空间 分 析 工 具 计 算 NDVI 和 气温 变量 的 
最 大 值 最 小 值 平均 值 和 标准 差 ,分析 其 时 间 变 化 特征 ,并 分 析 了 NDVIg 和 气温 的 变异 系数 (CV ,Coefficient 
of variation ) 。 
2.3 植被 变化 与 气温 、 热 通 量变 化 空间 交叉 验证 

为 了 量化 生长 季 (5 一 9 月 ) 黄 土 高 原 植被 生长 变化 对 气温 的 影响 ,本 文选 取 两 个 典型 时 段 进 行 差 值 分 析 : 
1998 一 2000 年 (黄土 高 原 退 耕 还 林 ( 草 ) 工程 初期 ,以 下 简称 S1) 与 2008 一 2010 年 (工程 开展 10 年 后 ,以 下 简 
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称 S2)。 首 和 完 计算 S1、S2 两 时 段 间 平均 年 NDVIg 气温 变量 及 热 通 量 的 变化 量 ,然后 通过 AreGIS 10.0 的 空间 
分 析 工 具 对 变化 量 进行 登 加 分 析 , 以 检测 植被 履 盖 度 与 气温 和 地 表 热 通 量 之 间 的 关系 ,并 探寻 其 影响 机 制 。 


3 结果 与 讨论 


3.1 植被 时 空 变化 

在 退耕 还 林 ( 草 ) 开 始 的 10 年 间 , 黄 土 高 原 NDVIg 整体 增加 ( 表 1) ,其 中 Sl 和 S2 期 时 段 NDVIg 均值 分 
别 为 0.38 和 0.45 ,增加 幅度 高 达 18.42% ;同期 ,NDVIg 最 大 值 也 增加 了 13.16% 。 而 且 ,S1 时 段 NDVIE 的 CV 
值 比 S2 高 ,表明 S1 时 段 区 域内 单元 的 NDVIg 差异 大 ,分 布 不 均 ,而 S2 时 段 的 区 域内 差异 较 小 。 


表 1 黄土 高 原 植 被 NDVIg 统计 值 
Table 1 The statistic characters of NDVIg on the Loess Plateau 


时 自 平均 人 最 大 人 最 小 值 标准 关 变异 系数 /% 
Time periods Mean Maximum Minimum Standard deviation Coefficient of variation 
Sl 0.38 0.83 0.02 0.17 44.74 
S2 0.45 0.88 0.01 0.18 40.00 


S1 :黄土 高 原 退 耕 还 林 ( 草 ) 工程 初期 ,the initial stage of the Grain for Green Project (GGP) ;S2; 工程 开展 10 年 后 ,10 years after GGP 


两 个 时 段 NDVIg 及 两 时 段 变 化 值 的 空间 分 布 差异 明显 (图 2) ,植被 恢复 区 域 明 显 增加 。 无 植被 覆盖 或 
植被 较 少 覆盖 区 域 (0<NDVIg 友 0.3) 从 S1 时 段 的 总 面积 占 比 39% 下 降 至 $2 时 段 的 25.10% ;植被 高 覆盖 区 域 
(NDVIg=0.6) 的 面积 从 5% 增 加 至 10.32%( 图 2a、 图 2b)。 和 黄土 高 原 中 部 \ 东 部 < 东南 部 的 绿地 覆盖 面积 增加 
较 快 ,NDVIg 增加 的 区 域 占 到 全 区 的 92.37% ;强烈 的 恢复 过 程 伴随 着 有 小 部 分 区 域 的 植被 退化 过 程 ( 图 2c)， 
在 兰州 以 北 较为 明显 ,NDVIg 下 降 最 低 值 达 到 -0.28。 退 化 区 域 可 能 受到 自然 状况 与 人 类 活动 的 共同 作用 ,该 
区 域 地 处 半 干 旱 偏 早 区 ,土壤 竣 薄 ,水 资源 短缺 ;生态 治理 难度 较 大 ” ;可 能 存在 的 过 度 放牧 也 会 对 干旱 区 风 
沙 草 原 植被 退化 有 着 重要 的 驱动 作用 。 


(a) S1 


(b) S2 
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图 2 黄土 高 原 不 同时 段 NDVIg 均值 及 两 时 上段 差 值 空间 分 布 
Fig.2 The spatial distribution of the mean NDVIg in 1998 一 2000 (S1) and 2008 一 2010 (S2), and the difference of NDVIg between S2 and 
Sl 
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3.2 ”站 点 尺度 植被 覆盖 与 气温 变化 

统计 得 到 54 个 站 点 S1、S2 时 段 NDVIg 与 气温 变量 统计 值 如 表 2 所 示 。 站 点 尺度 的 植被 覆盖 变化 与 像 元 
尺度 的 结果 一 致 ( 表 2 .图 2) ,NDVIg 各 统计 值 在 S2 时 段 均 出 现 上 升 ;同时 各 气温 变量 最 大 值 与 均值 下 降 。 其 
中 ,各 站 点 的 NDVIg 最 大 值 增加 幅度 最 大 , 达 0.07( 表 2) 。 研 究 时 段 内 ,生长 季 各 气温 变量 了 7T ，、 了 的 
均值 与 最 大 值 都 下 降 ,表明 降温 是 该 区 的 一 致 趋势 , 尤 以 有， 最 为 明显 。 与 之 相反 的 是 ,S1 至 S2 时 段 , 统 计 的 
气温 变量 的 最 小 值 都 升 高 。 这 说 明黄 土 高 原 气 温 状况 仍 受 到 大 太 度 全 球 变 暖 背景 的 影响 ,使 得 气温 变量 最 小 
值 出 现 上 升 ; 但 其 作用 不 如 局 地 强迫 明显 ,黄土 高 原 气 温 总 体 上 呈 下 降 趋 势 。 而 且 ,S1 时 段 NDVIg 与 各 气温 
变量 的 CV 值 比 S2 高 ,表明 S1 时 段 区 域内 单元 的 NDVI 差异 大 ,分 布 不 均 , 而 S2 时 段 的 区 域内 差异 较 小 。 进 
一 步 对 两 个 时 段 NDVI 以 及 气温 指标 的 均值 统计 量 进行 99% 置 信 区 间 的 配对 样本 TT 检验 ,NDVI、7,, Z 的 
均值 的 在 P=0.01 水 平 上 存在 显著 差异 ,了 ， 值 在 两 时 段 差 异 不 显著 。 


表 2 站 点 尺度 植被 NDVIg 及 气温 变量 统计 值 


Table 2 The NDVIg and temperature variable statistics in the Loess Plateau at site scale 


时 段 变量 单位 平均 值 最 大 值 最 小 值 标准 差 变异 系数 /% 
Time periods Variables Units Mean Maximum Minimum Standarddeviation Coefficient of variation 
S1 NDVIg / 0.42 0.75 0.12 0:19 45.24 
， Tt 19.16 23.98 9:31 3.04 15.87 
wj 25.61 30.29 15.44 .31 12.92 
全 13.46 18.87 4.23 3.03 22.51 
S2 NDVIg 0.49 0.82 0.13 0.20 40.82 
人 本 18.91 23.90 10.06 2.09 11.05 
a 于 25.28 29.89 16:26 3.11 12.30 
1 所 13.45 18.85 5.07 2.96 22.00 


A 
尘 
邮 


NDVIg: 植 被 生长 季 NDVI,NDVI in growing season;7, :气温 平均 值 ,Mean temperature;7 :气温 最 大 值 , Maximum temperature ;7 : 


小 值 , Minimum temperature 


研究 时 段 内 NDVIg 变化 (ANDVIg) 7 变化 ( AT ) 7 变化 (AT ss.) 与 7 变化 (ATu) 的 空间 分 布 如 
图 3 所 示 。 图 3 表明 气温 变化 的 空间 格局 与 NDVIg 变化 大 致 呈 相 反 模 式 ,植被 覆盖 增加 (减少 ) 的 站 点 气温 
下 降 ( 上 升 ) 。 在 退耕 还 林 ( 草 ) 开 始 的 10 和 年间, 黄土 高 原 中 部 偏 南 靠近 黄河 主干 道 的 站 点 NDVIg 增长 幅度 远 
大 于 其 他 地 区 站 点 ,最 大 的 在 离 石 站 ;增幅 达到 0.19。 研 究 时 段 内 黄土 高 原 植被 覆盖 整体 呈 增 加 趋势 ,少数 站 
点 的 NDVIg 出 现下 降 ,最 大 值 出 现在 黄土 高 原 东北 部 的 集 宁 站 ,平均 NDVIg 下 降 0.04。 黄 土 高 原 西北 部 站 点 
的 NDVIg 值 变化 较 小 。 

植被 覆盖 增加 地 区 的 气温 下 降 现象 比较 普遍 (图 3) ,各 站 点 了 值 在 中 部 偏 南 靠 近 黄 河 主干 道 下 降 情 况 
较 明 显 ,在 阳 果 站 达到 最 低 -1.1 % ;西南 部 有 少数 变 暖 现 象 , 变 暖 幅度 为 0 一 1 % ,最 高 的 增 温 出 现在 佛 坪 站 
0.79 % (网 3b)。A7， 值 分 布 与 AT， 值 大 致 相同 ,西南 部 升温 中 部 降温 ,分 别 在 玛 曲 站 (0.82 %C ) 与 榆林 站 
(一 1.95-%C ) 达 到 最 高 (图 3c)。AT7T,, 空间 分 布 模式 与 AT 与 AT 分布 模式 存在 差异 ,其 在 中 部 地 区 有 少数 
变 暖 站 点 ;幅度 在 0 一 0.5%C , 变 暖 站 点 仍 集中 分 布 于 西南 地 区 ,最 高 的 增 温 出 现在 佛 坪 站 达 1.25% ;降温 最 低 
的 站 点 西 移 至 固原 北部 的 盐池 站 ,达到 -1.54 % (图 3d)。 黄 土 高 原 及 周边 植被 覆盖 站 点 的 日 最 低 气 温 空 间 
分 布 变 异性 高 于 日 平均 气温 与 日 最 高 气温 。 

对 研究 时 段 内 平均 NDVIg 气温 变化 的 交叉 验证 分 析 表 明 ,植被 覆盖 增加 与 气温 变量 降低 在 空间 上 呈正 
相关 ,上 且 植 被 对 气温 的 影响 在 日 间 较 夜间 大 (图 3、 图 4)。 在 研究 时 段 内 7, 7 的 第 区 象限 站 点 数 最 多 , 占 
黄土 高 原 的 大 部 分 区 域 ,分 别 约 占 64.81%、72.22%;Ti, 第 1 象限 站 点 数 大 致 相同 ,分 别 约 占 44.44%、 
46.30%。 这 表明 植被 对 气温 的 影响 在 日 间 较 夜间 大 ,可 能 是 因为 植被 与 大 气 间 的 大 部 分 能 量 交 换 活 动 都 在 
有 光照 的 日 间 。 站 点 数 在 第 亚 象限 空 间 上 呈 极 少数 分 布 ,了 .7 .与 Ti 站 点 数 分 别 占 5.56%、3.70% 和 
5.56% 。 
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图 3 基于 气象 站 的 NDVIg、ThensTuu 与 Tim 生长 季 变 化 情况 


Fig.3 Site-based Changes of NDVIg, mean temperature (Ten ), maximum temperature ( 7、 ) ，and minimum temperature (Tm) in 


the growing season ( mean of 1998 一 2000 minus mean of 2008 一 2010 ) 


-3.3 ”地表 热量 平衡 分 析 
< 上 述 结果 表明 ,植被 覆盖 增加 可 能 会 对 黄 士 高 原 的 局 地 热量 状况 造成 影响 , 带 来 一 定 的 降温 效应 。 植 被 
可 以 通过 改变 下 垫 面 的 反照 率 .粗糙 度 及 蒸 散 等 方式 影响 大 气 ,从 而 影响 地 面 环流 和 能 量 平 衔 , 引 起 局 地 气温 
l 的 变化 一 。 为 进一步 验证 研究 结果 ;使 用 热量 平衡 分 析 的 方法 对 植被 覆盖 变化 对 气温 影响 的 基本 过 程 及 其 
-= ”可 能 的 物理 机 制作 一 些 探讨 ,选择 影响 地 表 能 量 平衡 关系 的 物理 量 '” :地 - 气 之 间 的 感 热 (ASH) 、 潜 热 通 量 
二 (ALH) 与 大 气 下 行 长 波 辐射 (ALD) -另外 ,还 对 研究 时 段 内 生长 季 植 被 覆盖 度 变化 分 级 为 3 类 :高 绿化 站 点 
(ANDVIg>0.1) .中 绿化 站 点 (0.05<ANDVIg 大 0.1) 与 低 绿化 站 点 (0<ANDVIg 大 0.03 ) ,定量 统计 其 对 应 气温 与 
热 通 量变 化 值 如 表 3 所 示 

各 气温 变化 与 ANDVIlg 变化 高 度 一 致 ( 表 3) ,气温 变化 对 应 NDVIg 绿化 区 间 呈 梯度 分 布 , 即 随 着 绿化 程 
度 的 增加 ,各 气温 变量 的 降温 程度 也 相应 增加 。7,, 变 量 在 中 绿化 与 低 绿 化 站 点 出 现 少 量 增 温 。 除 ALD 在 低 
绿化 站 点 下 降幅 度 较 中 绿化 站 点 高 外 , 热 通 量变 化 与 ANDVIg 变化 相 一 致 , 随 着 绿化 程度 增加 ,ALH 增加 的 幅 
度 增加 ,ASH 下 降 的 幅度 增加 ,ALD 下 降 的 幅度 增加 。 


表 3 S1 至 S2 时 段 不 同 ANDVIg 增长 幅度 对 应 气温 变量 与 热 通 量变 化 
Table 3 Mean temperature variables and heat fluxes changes at sites grouped by the amplitude of NDVIg increased ( ANDVI) between S1 


and S2 
区 间 Range ANDVIg AT LSC AT /CC AT ZL%C a pe i 
>0.1 0.13 -0.44 -0.10 -0.64 14.92 -13.52 -5.72 
0.05 一 0.1 0.07 -0.26 0.02 -0.37 6.70 —4.44 -5.14 
0 一 0.05 0.03 -0.11 0.05 -0.14 3.98 -3.75 -6.15 


ANDVIg: 植 被 生长 季 NDVI 变化 ,NDVI change in growing season;AT :人 气温 平均 值 变化 , Mean temperature change;7,, :气温 最 大 值 变 化 ， 
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ANDVIg 与 Tinean 的 关系 
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NDYVI change in growing season 


图 4 黄土 高 原 植被 NDVJg 与 气温 变化 关系 
Fig.4 The relationships between the changes ofthe temperature and NDVIg on the Loess Plateau 


Maximum temperature change;7 :气温 最 小 值 变 化 ,Minimum temperature change; ALH : 潜 热 通 量变 化 ,Latent heat flux change; ASH: 感 热 通 量变 化 ， 


Sensible heat fluxchange; ALD :下 行 长 波 辐射 变化 ,Downward longwave radiation cliange 


为 探讨 ANDVIg 与 气温 变化 之 间 的 物理 关联 ;交叉 验证 了 热 通 量变 化 与 NDVIg 变化 之 间 的 关系 。 由 图 5 
可 得 ,在 退耕 还 林 ( 草 ) 开 始 的 10 年 间 , 黄 主 高原 植 被 攻关 增 加 与 洲 热 通 量 增加 、 感 热 通 量 与 下 行 长 波 辐 射 下 
降 在 空间 上 呈正 相关 关系 ( 表 3, 网 $)z 植被 覆 善 增加 后 ,地 表 潜 热 通 量 由 于 区 域 土 壤 含水 量 的 增加 、 地 表 粗 
糙 率 及 叶 面 积 指 数 的 增 大 及 节 散 发 的 增加 而 出 现 了 显著 增加 , 感 热 通 量 减少 。 如 图 5 所 示 , 约 有 64.81% 的 站 
点 植被 履 盖 与 潜 热 通 量 同时 增加 ,增加 幅度 为 0.46 一 44.77 W/m ; 约 为 74.07% 的 站 点 植被 覆盖 增加 但 感 热 通 
量 减少 ,减少 幅度 为 =0.22 一 30.53 W/m 。 

长 波 辐射 变化 可 以 从 植被 地 面 反 照 率 方面 进行 解释 ,一 般 而 言 ,植被 的 反照 率 比 荒漠 地 要 低 得 多 。 黄 
土 高原 地 区 植被 覆盖 增加 后 引起 地 表 反 照 率 的 改变 ,会 对 地 表 吸 收 短波 太阳 辐射 和 发 射 长 波 辐射 造成 影响 。 
许多 人 研究 证 明 ;植被 覆盖 的 增加 ,可 以 通过 减少 地 表 反 照 率 ,增加 地 表面 吸收 的 太阳 辐射 量 ,从 而 产生 对 气温 
的 正 反馈 然而 ,这 种 正 反 照 率 反馈 可 能 部 分 被 更 高 的 蒸 散 发 的 冷却 作用 所 抵消 。 这 种 蒸 散发 带 来 大 
气 中 水 汽 增 加 ,使 得 云 量 增加 ,会 造成 下 行 长 波 辐射 的 明显 减少 ,如 图 5 所 示 ,减少 的 站 点 分 别 占 黄土 高 原 的 
79.63% ,幅度 在 -0.02 一 13.91 W/m?。 

从 热量 平衡 的 角度 分 析 ,植被 增加 使 得 土壤 含水 量 增加 ,表面 阻抗 增加 ,植被 的 蒸腾 增加 ,以 及 粗糙 度 增 
加 而 造成 的 泊 流 垂直 扩散 增强 ,会 造成 潜 热 通 量 增加 , 感 热 通 量 的 减少 。 感 热 通 量 减 小 幅度 小 于 潜 热 的 增加 
幅度 , 故 该 区 的 有 效 热 通 量 ( 感 热 + 感 热 ) 是 增加 的 ,数值 在 0.24 一 14.24 W/m 。 再 加 上 反照 率 减少 ,吸收 的 净 
短波 辐射 增加 。 通 过 分 析 热 量 平 衡 分 量变 化 ,由 于 地 表 净 吸收 的 辐射 增加 , 且 有 效 热 通 量 的 值 增加 ,要 在 更 低 
的 气温 下 才能 使 地 表 吸 收 的 短波 、 发 射 的 长 波 、 输 送 的 有 效 热 通 量 达 到 新 的 平衡 。 通 过 植被 燕 散 作 用 加 强 并 
可 能 导致 降温 的 结果 与 梁 玲 等 ” 、 苟 娇 娇 等 ”在 黄土 高 原 进 行 的 植被 变化 对 气候 影响 试验 数值 模拟 研究 一 
致 。 这 与 北极 地 区 植被 增长 反馈 作用 分 析 ' 及 中 国 的 土地 覆盖 变化 敏感 性 试验 研究 结果 不同, 他 们 认为 
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ANDVIg 与 ADLH 的 关系 
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图 5 黄土 高 原 植被 NDVIg 变化 与 热 通 量变 化 关系 
Fig.S The relationship between the changes'in the heat fluxes and NDVIg on the Loess Plateau 


植被 增长 使 得 地 表 反 照 率 减 少 , 增 加 了 表面 吸收 的 太阳 辆 射 量 ,进而 导致 植被 对 气温 的 正 反 饿 。 然 而 ,从 局 地 
矿 度 上 来 看 ,黄土 高 原 气 候 冬季 严寒 .夏季 暖 热 , 且 大 量 履 盖 草 地 罗 被 。 植 被 的 蒸 散发 作用 在 该 地 区 有 者 重要 
的 地 位 ,其 可 以 在 植被 增加 的 情况 下 削弱 由 于 反照 率 减 少 导 致 的 变 暖 ,植被 需要 蒸腾 更 多 的 水 分 ,从 而 在 生长 
季 持 续 降温 反 饿 。 植 被 恢复 可 能 会 对 黄 士 高原 气候 变 暖 现象 有 减缓 作用 。 


4 结论 


本 文 基于 气象 观测 数据 与 卫星 遥感 NDVI 数据 ,对 黄土 高 原生 态 工程 实施 前 后 10 年 间 (S1:1998 一 2000 
年 ,工程 初期 ;$2 :2008 一 2010 年 ,工程 T0 年 后 ) 的 植被 增长 状况 与 气温 变量 时 空 变 化 状况 进行 交叉 验证 , 探 
讨 黄土 高 原 植 被 恢复 对 气温 的 反馈 作用 ;同时 基于 EIN-Interim 再 分 析 资 料 提取 的 地 面 观测 热 通 量 值 变化 ,分 
析 植 被 变化 对 气温 影响 的 基本 过 程 及 其 可 能 的 物理 机 制 , 探 讨 植 被 与 大 气 间 相互 作用 。 通 过 研究 得 出 以 下 
结论 : 

(1) 黄 卡 高 原 植被 恢复 比较 明显 ,生长 季 植 被 增加 显著 ; 与 之 对 应 的 是 生长 季 3 个 气温 变量 (了 7 、 
7,, ) 的 均值 与 最 大 值 在 S2 时 有 段 都 下 降 , 尤 以 7, 最 为 明显 。 

(2) 植 被 覆盖 增加 与 气温 变量 降低 在 空间 上 呈正 相关 。 同 时 ,植被 覆盖 增加 与 潜 热 通 量 增加 、 感 热 通 量 
与 下 行 长 波 辐射 下 降 在 空间 上 也 呈正 相关 关系 ,植被 恢复 通过 加 剧 植被 蒸 散 作用 导致 近 地 面 降温 。 

研究 还 存在 一 定 不 足 ,数据 的 不 确定 性 和 模型 缺陷 会 阻碍 对 植被 活动 增加 引起 的 气温 反馈 大 小 的 定量 结 
论 寺 未 来 应 考虑 加 大 观测 数据 的 测量 时 空 精度 .在 模型 模拟 中 引入 其 他 参数 (如 叶 面 积 指数 .表面 反照 率 ) 使 
陆 面 过 程 参数 化 更 符合 实际 状况 。 这 些 改进 预期 能 够 更 好 地 了 解 植 被 在 气候 系统 中 的 作用 。 
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